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Abstract 

The review contains an analysis of the observed and model curves of the 
interstellar extinction and polarization. The observations mainly give in- 
formation on dust in diffuse and translucent interstellar clouds. The fea- 
tures of various dust grain models including spherical/non-spherical, homo- 
geneous/inhomogeneous particles are discussed. A special attention is de- 
voted to the analysis of the grain size distributions, alignment mechanisms 
and magnetic field structure in interstellar clouds. It is concluded that the 
interpretation of interstellar extinction and polarization is not yet complete. 

Keywords: Light scattering, Nonspherical particles. Composite particles. 
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1. Introduction 

The properties of cosmic dust grains in various objects from comets to 
distant galaxies are derived from observations of interstellar extinction, in- 
terstellar polarization, scattered radiation, infrared (IR) continuum emission 
and IR features. Modelling of these observations is aimed at estimates of 
the grain size, chemical composition, shape, structure, and alignment. The 
observed wavelength dependencies of interstellar extinction and polarization 
(interstellar extinction A{X) and polarization P{X) curves) still remain the 
main sources of information on dust in diffuse and translucent interstellar 
clouds. 



Email address: nvvSastro . spbu . ru (N.V. Voshchinnikov) 



Preprint submitted to JQSRT 



March 8, 2013 



In this review, we discuss observations and modelling of the interstellar 
extinction and polarization curves. Special attention is paid to various dust 
grain models. For an extended consideration of the properties of dust in 
different astronomical objects see [H, 0, S, 0, [ll- Excellent historical reviews 
on dust astrophysics are given by Dorschner ]6[ and Li . 



2. Extinction 

2.1. Extinction curve: production and fitting 

The normalized extinction curves A^'^\\^^) commonly studied can be 
calculated as the ratio of the colour excess of a star behind the dusty cloud 
E{X — V) to the colour excess E{B — V) 

Here, As and Ay is the extinction in the B and V bands, respectively. Some- 
times, the normalization of the extinction curve on Ay or extinction in an- 
other band is used. 

In such a manner we can determine only the "selective" extinction (red- 
dening), i.e. the difference of extinction at two wavelengths. The absolute 
value of extinction can be found as 

Av = RvE{B-V), (2) 

where the coefficient Ry is often evaluated from observations in the visible 
and IR taking into account that Ax ^ when A — t- oo. From Eq. ([T]), it 
follows that 

R - - - ^) (o^ 

- E{B -V)~ E{B - V) • 



In the diffuse medium on average Ry = 2.4 — 3.6 [8||. 

In general, the interstellar extinction curve has a power law-like rise from 
the IR to the visible, a prominent feature (bump) near A 2175 A, and a steep 
rise in the far-UV. This rise is a manifestation of the very strong feature 
with a maximum near A ~ 700 A j9|. Figure P shows the extinction curve 
averaged over 243 Galactic B and late-0 stars |8|. In the near IR- visible part 
of the spectrum, the distinction between the extinction curves of different 
stars is rather small. The IR extinction at wavelengths A = 0.7-5 fim was 
approximated by the power-law dependence: A{X) oc A"'' with /3 = —1.84 
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Figure 1: The average extinction curve for 243 Galactic stars with 2.4 < Ry < 3.6 and the 
extinction curves in the direction of six stars in Sco-Oph Q. The values of the coefficient 
Rv and the UV fitting coefficients C4 and C5 (Eq. ([4])) are indicated in the legend. The 
effect of variations of coefficients C4 (left panel) and C5 (right panel) is illustrated. Adapted 
from 81. 



Fitzpatrick and Massa 11 1 found that the index (3 was i?y-dependent: 
/9 > 2 if Rv < 3 and (3 < 1.5 if Ry > 3. The mid-IR extinction at wavelengths 
A ^ 3 fim measured in the Galactic plane 12| and toward the Galactic centre 



13l | becomes grayer than the near-IR extinction. 

In the UV region, the extinction curves differ strongly [s], Ei, lif (see also 



Fig. H]) demonstrating that the mean curve in the UV obviously has little 
meaning. The position of the UV bump center varies a little from star to 
star and occurs at Aq = 2174 ± 17A or Aq^ = 4.599 ± 0.012 /im'^ The 
total half-width of the bump is = 0.992 ± 0.058 /xm"^ that corresponds to 



470±27A [16 

Dorschner [17[ was the first who suggested to approximate the shape 
of the bump profile by the classical (Lorentzian) dispersion profile. After 
examination of the lUE extinction curves for many lines of sight, Fitzpatrick 
and Massa js], [igI, 18 | deduced a single analytical expression with a small 
number of parameters describing extinction in the region 1150 A < A < 
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2700 A {x = X-^) 



X] 



E{X~V) 
E{B-V) 



ci + C2X + csD{x, W, xo), 

ci + C2X + csD{x, W, xo) + ca{x - 05)^, 



X < C5, 
X > C5, 



where 



D{x,W,Xo) 



X 



[X^ 



)2 +xW2' 



(4) 
(5) 



Originally, the value of C5 was fixed at 5.9 /xm~^ [18|. Equation consists of: 
i) a Lorentzian-like bump term (requiring three parameters, corresponding 
to the bump width W, position xq, and strength C3), ii) a far-UV curvature 
term (two parameters C4 and C5; see Fig. [T]for illustration), and in) a linear 
term underlying the bump and the far-UV (two parameters ci and C2). The 
parameters of the average extinction curve presented in Fig. [1] are: Ry = 
3.001, xo = 4.592 /im-i, W = 0.922 /im"!, ci = -0.175, C2 = 0.807, C3 = 
2.991, C4 = 0.319, C5 = 6.097 jsj. 

Fitzpatrick and Massa [sj] note that there is no correlation between the 
UV and IR portions of the Galactic extinction curves. This fact is illustrated 
by Fig. |2] where the dependence of the coefficient Ry on the quantity describ- 
ing the strength of the far-UV extinction curvature A1250 = 04(8.0 — cs)^ is 
shown. Absence of the correlation between Ry and A1250 contradicts the 
often used representation of the extinction curves from the UV to IR as a 
one-parameter family dependent on Ry (so called CCM model introduced 
by Cardelli, Clayton and Mathis 19,[2o[). As mentioned in jsl, the relations 



between Ry and UV extinction found in [19|, |20| can arise from sample selec- 
tion and methodology. Evidently, some bias relates to the number of clouds 
available in the line of sight, i.e. to the distance to the star. There exists a 
wide scatter of the data for nearby stars in comparison with the distant stars 
(see Fig. A large difference in extinction for the nearby stars observed 



through single clouds was first noted by Krelowski and Wegner [2l| . It should 



be mentioned that the major part of anomalous or peculiar extinction sight- 



lines studied so far are related to not very distant stars |22l . l23l |24J . l25l |26 



2.2. Interpretation: homogeneous spheres 

The extinction of stellar radiation at the wavelength A after passing a 
dust cloud is equal to 



A{X) = -2.5 log /(A)//o(A) ^ 1.086r(A), 



(6) 
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Figure 2: The coefficient Ry in dependence on the quantity A1250 = 04(8.0 — 05)^ showing 
the strength of the far-UV extinction curvature in the direction of 321 stars with known 
distances from 8] . Fihed and open circles show data for stars with the distances D < 1 kpc 
and D > 1 kpc, respectively. Cross corresponds to the average Galactic extinction curve. 

where Io{X) is the source (star) intensity, r(A) the optical thickness which 
can be found as the total extinction cross-section of all particles types along 
the line of sight in a given direction. Interpretation of the interstellar ex- 
tinction is often performed using homogeneous spherical particles of various 
size distributions. Then, the wavelength dependence of extinction can be 
calculated as 



A{X) = 1.086 ^ ^ j Cext,i(m, r„ A) n,(r,) dr,j dl. 



(7) 



Here, nj{rs) is the size distribution of spherical dust grains of the type 
j and radius with the lower cut-off rs^min and the upper cut-off rg max, 
Cext(---) = Trr^ (5ext(---) the extinction cross-section, Qext the extinction efii- 
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ciency factor, m the refractive index, D the distance to a star. From Eq. ([7]), 
an important conclusion follows: the wavelength dependence of interstellar 
extinction is completely determined by the wavelength dependence of the 
extinction efficiency Qext- 




Figure 3: Wavelength dependence of the extinction efficiency factors for homogeneous 
spherical particles of different sizes consisting of astronomical silicate and amorphous car- 
bon. The dashed segment shows the approximate wavelength dependence of the mean 
Galactic extinction curve at optical wavelengths. Adapted from [2]. 

The average interstellar extinction curve in the visible-near UV can be 
approximated by the power law A{X) oc A~^ '^ [8]. Such wavelength depen- 
dence can be produced by submicron-sized particles of the typical radius 
(r) ^ 0.05 — 0.1 /xm. In this case, for more absorbing materials like amor- 
phous carbon or iron we have smaller particles and for less absorbing ma- 
terials like silicate or ice we need larger particles (see Fig. |3]). So, from the 
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wavelength dependence of extinction only the product of the typical particle 
size and the refractive index (r) \m — 1\ ^ const, can be determined, but 
not the size or chemical composition of dust grains separately. In order to 
solve this problem, the dust-phase abundances of the main elements forming 
dust (C, O, Mg, Si and Fe) need to be taken into account and to reproduce 
the absolute extinction. Unfortunately, despite numerous observations of the 
interstellar absorption lines (see a compilation of Gudennavar et al. (j?! ) 
abundances with good accuracy are known just for a restricted number of 
diffuse and translucent clouds |28|, |29|, |30 . 

Another problem having many solutions is identification of the UV bump 
near A 2175 A. Various materials with isotropic and anisotropic properties 
such as silicate (enstatite), irradiated quartz, oxides (MgO, CaO), organic 
molecules have been considered as carrier candidates (see discussion in [l|,l2|). 
However, the position and width of the bump are strongly suggestive of 
TT — 7- TT* transitions in graphitic or aromatic carbonaceous species domi- 
nating by sp^ bonding. Therefore, small graphite particles and polycyclic 
aromatic hydrocarbons (PAH molecules) are considered as the favourite ma- 



terials [311 |32|, |33|. Unfortunately, the reliable identification of the carrier 



remains unknown. The attempts to find the UV or visual bands of PAHs 



have failed |34j . |35| . We cannot even determine the size of graphite spheres 
responsible for the A 2175 A feature (see Fig.H]). The profiles with the central 
position near Aq ^ = 4.6 /im~^ can be obtained if we take particles with the 
radius ^ 0.015 /im (Fig. IH left panel). Although for single-size particles 
the width of the calculated profiles is smaller than the observed one, a simple 
bi-modal size distribution allows a fit to both the position and the width of 
the mean Galactic profile (Fig. |U right panel). 

The far-UV extinction can be explained by tiny particles of the typical 
radius (r) ^ 0.01 — 0.03 fim (see Fig. [3]). The number density of such grains 
is ~ 1000 times large than the submicron particles producing the visual- 



near-IR extinction 36 



Because of temperature fluctuations, such particles 
are protected from growth by accretion in the interstellar clouds. The far- 
UV rise of extinction may be also fitted as the low-energy side of cr — )■ a* 



transitions in PAHs (see 37(| for discussion). 



By using particles of different chemical composition and applying Eq. ([7]) 
it is possible to interpret the interstellar extinction and to reconstruct the 
dust size distribution. In the pioneer work of Oort and van de Hulst 38 
the size distribution of icy grains was found in the tabular form. Later, 
Greenberg 39| fitted it with an exponential function. By using minimiza- 
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Figure 4: Normalized extinction efficiencies for graphite spheres. The curve marked as 
"observations" corresponds to the wavelength dependence of the UV bump given by the 
mean Galactic extinction curve. The central position of the observed UV bump and its 
range of variations are marked. The left panel shows extinction of single size graphite 
spheres. The right panel shows the summary extinction of two graphite spheres with 
radii Vs — 0.005 /im and rg = 0.03 /xm (from left panel) taken in equal proportions. All 
calculations were made in the "2/3-1/3" approximation for the averaged extinction factors 
Qoxt — 2/3Qoxt(e±) + 1/3 Qcxt(e||), where e± and e|| are the dielectric functions for two 
cases of orientation of the electric field relative to the basal plane of graphite. Adapted 
from [3]. 



tion of the statistic, Mathis et al. j4l| reconstructed the power-law size 
distribution for sihcate and graphite particles. These two simplest size distri- 
butions contain the only parameter (except for the lower and upper cut-offs, 
see Tabled]). More complicated two-parameter distributions were applied by 
Wickramasinghe and Guillaume 42[ and Wickramasinghe and Nandy 43 1 in 
order to fit the mean extinction curve with graphite grains and a mixture of 
graphite, iron and silicate grains, respectively. A comprehensive discussion 
of the early attempts to model the interstellar extinction can be found in the 
review of Wickramasinghe and Nandy fi^l . 

The size distributions of dust grains can be also found from extinction 
measurements by solving the inverse problem. It has been made by the 
maximum en trop y method 
tion 



45 



47, 48, 49 



and by Tikhonov's method of regulariza- 
The obtained size distributions were approximated by 
rather complicated functions containing up to 14 p arameters (see Tabled]). 
Complexity of the size distributions used in 3l|, |47| is explained by the chal- 
lenge in reproducing the diffuse Galactic IR emission as well. 
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Table 1: Dust size distributions used for interpretation of interstellar extinction. 
Author(s) (year) reference; size distribution function A^parametcrs 
Greenberg (1968) [39]; exponential 

n(rs) oc exp [-5 (rjr^of] 1 
Isobe (1973) [40]; exponential 

n(rs) cx exp [- (rs/r^o)] 1 
Mathis et al. (1977) [41]; power-law; MRN mixture 



n(rs) oc 



Wickramasinghe and Guillaume (1965) [42]; normal 
n(rs) oc exp [- {r^ - f^f /{2a^)] 

Wickramasinghe and Nandy (1971) [43]; lognormal 
n{rs) oc exp [-1/2 {rjr2) ] 

Kim et al. (1994) [45]; power-law with exponential decay 
n(rs) oc r~'^ exp(-rs/rsb) 

Mathis (1996) [46]; power-law with exponential decay 
n(rs) oc r~'^° exp[-(7irs + 72/rs + 'Js/r^)] 

Weingartner and Draine (2001) [31]; two lognormal 

n(rs) oc Vcir^) + " ' ' 



rs vn;C,Si 
l + /3c,Sirs/rt;c,Si, /? > 
(l-/3c,Sirs/rt;C,Si)-S /9 < 



exp{-[(rs -rt;c,si)At;C,Si]^}, > n-csi 
Zubko et al. (2004) [47]; 
logn(rs) = Co + 6olog(rs) - 6i| log(rs/ai)p^ 

-62I log(r,/a2)r2 - - a^r' " &4|rs - a,r^ 14 
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2.3. Interpretation: inhomogeneous and composite particles 

Progress in observations, the light scattering and grain growth theories 
gave rise to new dust models with grains more complicated than homogeneous 
spheres. 

Wickramasinghe 50|, l5l| was the first who studied the optical properties 
of core-mantle grains which could grow in interstellar clouds due to accre- 
tion of volatile elements on refractory particles. He calculated extinction 
produced by graphite core-ice mantle and silicate core-ice mantle spheres. 
Extensive calculations of extinction for graphite core-ice mantle particles 



were also made by Greenberg j39| who later proposed the existence of parti- 
cles with silicate cores coated by a layer of organic material in diffuse clouds 
and silicate-organic-ice grains in molecular clouds 36|. Such grains were a 



component of the dust mixture reproducing interstellar extinction [52 



The growth of interstellar grains due to their coagulation in dense molec- 
ular cloud cores may result in formation of grain aggregates with large voids 
53|. The internal structure of such composite grains can be very compli- 



cated, and their optical properties cannot be described by the model of core- 
mantle spheres. Exact calculations are possible for complex aggregates of 



rather small sizes [5J, |55|, |56|. Therefore, very complicated particles are re- 
placed by more simple "optically equivalent" ones. A very popular approach 
is to make calculations using the Mie theory for homogeneous spheres with an 
average refractive index derived from one of the mixing rules of the effective 
medium theory (EMT; see, e.g., |46|, |47|, |57| and Table |2]). Another possibil- 
ity to treat composite aggregate grains is to consider multi-layered particles. 
As shown by Voshchinnikov and Mathis 
and Voshchinnikov 



63 



60| for spheres and by Farafonov 



for spheroids, the scattering characteristics of layered 
particles slightly depend on the order of materials and become close to some 
"average" ones, when the number of layers exceeds 15 — 20. According to 
estimates made in the optical properties of layered particles resemble 
those of heterogeneous particles having inclusions of various sizes while the 
EMT-Mie approach can be used if the particles have small (in comparison 
with the wavelength of the incident radiation) "Rayleigh" inclusions. 

Inhomogeneous and composite particles have an advantage over homo- 
geneous ones as there exists the possibility of including vacuum as one of 
the materials. The new dust models with fluffy, porous particles are able to 
produce the same extinction with a smaller amount of solid material than 
dust models with compact particles. The amount of vacuum in a particle 
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Table 2: Models of inhomogeneous spherical grains used for interpretation of interstellar 
extinction. 



Author (s) (year) reference 

Wickramasinghe (1963) [50] : 
Greenberg (1968) [39] 

Wickramasinghe (1970) [51] 

Greenberg, Li (1996) [52] 

Mathis and Whiffen (1989) [57]; 
Mathis (1996) [46] 

Zubko et al (1998, 2004) [47, 49] 

Vaidya et al. (2001) [58] 

Voshchinnikov and Mathis (1999) [60]; 
Voshchinnikov et al. (2006) [6]J 

latiet al. (2008) [62]; 
Cecchi-Pestellini et al. (2010) [37]; 
Zonca et al. (2011) [24] 

Rai and Rastogi (2012) [26] 



Model 

graphite core-ice mantle 

silicate core-ice mantle 

silicate core-organic mantle 

EMT-Mie: silicate + amorphous 
carbon + iron + voids 

EMT-Mie: silicate + organic 
refractory + water ice + voids 

silicates with graphite inclusions 

multi-layered: vacuum/silicate/ 
amorpous carbon 

four- layered: vacuum- 

silicate-sp^-carbon- 

sp^-carbon 

nanodiamonds coated by 
amorpous carbon or graphite 



can be characterized by its porosity V {0 <V < 1) 

'P = Kac/Kotal = 1 — Kiolid/Vtotal- (§) 

The role of porosity in extinction is seen from Fig. |5]that gives the wavelength 
dependence of the normalized cross section 

Cext (porous grain) 

Cext (compact grain of same mass) 

/3 Qcxt (porous grain) 

<5cxt (compact grain of same mass) 

This quantity shows how porosity influences the extinction cross section. As 
follows from Fig. |5l as P increases the model predicts a growth of extinction 
of porous particles in the far-UV and a decrease in the visual-near-UV. In 



a 



(n) 



[i-rr 
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Figure 5: Wavelength dependence of the normahzed extinction cross section (see Eq. ^) 
for multi-layered spherical particles with compact = 0.1 /xm. The particles are of the 
same mass but of different porosity. If Cg"^ > !> extinction of a porous particle is larger 



than that of a compact one of the same mass. Adapted from 61 [ 



comparison with compact grains, layered particles can also produce rather 
large extinction in the near-IR. This is especially important for the expla- 
nation of the flat extinction across the 3 — 8 fim wavelength range measured 



for several lines of sight (see |l3|, |6l| for discussion). It is also seen from 



Fig. |5] that an addition of vacuum into particles does not lead to a growth 
of extinction at all wavelengths and material saving. Evidently, the final 
conclusion can be made after detailed comparison of the observations with 
theoretical calculations at many wavelengths. 

Table |2] contains information about models of inhomogeneous spherical 
grains used for the interpretation of interstellar extinction. Inhomogeneous 
non-spherical particles of the simplest shapes (cylinders, spheroids) are also 
used for simultaneous interpretation of interstellar extinction and polariza- 
tion (see Table H]). In this case major attention has been paid to the modelling 



^This supports a conclusion of Li [65[ that an interpretation of the observed interstellar 
extinction curve using only very porous grains should not give any gain in dust-phase 
abundances. 
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of polarization because extinction usually has only a slight dependence on 
the particle shape and orientation (2I. lod^. 

The models discussed in this Section were first applied to interpretat- 
ing the average Galactic extinction curve. Modelling of extinction in several 



sighthnes has been also performed |22|, |23i |26|, |3l|, |48|, |49|, |6ll |67|, |68|, l69|, [70 



The models include multi-component mixtures of bare particles with a rather 
complicated size distribution function or/and inhomogeneous particles and 
PAHs and often took into account abundance restrictions. Especially popu- 
lar is the direction to the halo star HD 210121 with very high UV extinction. 
The extinction for this star has been modeled by Li and Greenberg jloj with 
a mixture of silicate core-organic mantle spheres, bare spheres and PAHs, by 
Larson et al. 70| and by Clayton et al. |7l| with two-component (silicate, 
graphite) and three-component (silicate, graphite, amorphous carbon) mod- 
els, by Weingartner and Draine 3l| who considered the mixture of carbona- 
ceous and silicate grains with the size distribution given in Table [H and by 
Rai and Rastogi 26| who used a silicate-graphite mixture and nanodiamonds 
coated by carbon. This list illustrates the non-uniquiness of parameters of 
dust grains obtained from modelling. 

A similar conclusion about the ambiguity of the modelling follows from 
several attempts to interpret the peculiar extinction curves characterized by 
a broad A 2175 A bump and a steep far-UV rise or a sharp A 2175 A bump 
and fiat far-UV extinction (see Fig. [T]). Using the model of Weingartner and 
Draine j3l| (see Tabled]), Mazzei and Barbaro 23 1 derived the parameters of 
the size distribution for 64 stars. Variations of the parameters were attributed 
to the selective grain destruction in both shocks and grain-grain collisions. 



Zonca et al. [2J] found an excellent fit for different extinction curves for 
15 sightlines with the mixture of layered porous grains for reproduction of 
the near-IR- visual extinction and PAHs to account for the bump and far-UV 
extinction. Large grains consisting of silicates coated by layers from graphitic 
and polymeric amorphous carbons (see Table [2]) were suggested in the model 
of Jones et al. 72[ (see 73 1 for more details). Rai and Rastogi 26 1 analyzed 
anomalous extinction curves in the direction of 10 stars and showed that a 
very good match with the far-UV rise of extinction was obtained if to include 
nanodiamonds coated by graphite or amorphous carbon as a component of 
the silicate-graphite mixture. 

Summarizing this discussion of the interstellar extinction, it should be 
noted that there is a wide diversity in the models and the non-uniqueness in 
the results. In spite of numerous attempts to use very complicated inhomo- 
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geneous particles, the Mie theory for homogeneous spheres keeps its leading 
position as a main tool for the interpretation of interstellar extinction. Fur- 
ther progress in the investigations should include a clear role for PAHs and 
modelling of the extinction on the basis of interstellar abundances in selected 
directions. A sofisticated undertstanding of the origin of UV extinction, the 
available solid-state material and the grain growth process would stimilate 
going from simple Mie theory to justified models of complex particles. 



3. Polarization 



3.1. Observations: Serkowski curve and polarizing efficiency 

Interstellar linear polarization is caused by the linear dichroism of the in- 
terstellar medium due to the presence of non-spherical aligned grains. Dust 
grains must have sizes close to the wavelength of the incident radiation and 
specific magnetic properties to efficiently interact with the interstellar mag- 
netic field. The direction of alignment must not coincide with the line of sight 
and there must be no cancellation of polarization during the propagation of 
radiation through the interstellar medium. 

Interstellar polarization was discovered in 1949 by Hiltner 



74|, Hall [75 



and Dombrovskii 76|| in the course of the search for the Sobolev-Chandrase- 
khar effectl^. The wavelength dependence of polarization P(A) in the vis- 
ible part of spectrum is described by an empirical formula suggested by 



Serkowski 77 



P(A)/Pmax = exp[-inn'(A^ax/A)]. 



(10) 



This formula has three parameters: the maximum degree of polarization 
Pmax, the wavelength corresponding to it Amax and the coefficient K char- 
acterizing the width of the Serkowski curve. Initially, the coefficient K was 
chosen by Serkowski 77, 78| to be equal to 1.15*^1^. 



^Sobolev and Chandrasekhar have shown that the polarization of radiation at the hmb 
of a star due to the electronic (Thomson) scattering should reach ^ 12%. Eclipsing binaries 
with extended atmospheres have been suggested in order to observe the effect. 

■^The Serkowski curve is just one of possible approximations of the observed dependence 
P(A). For example, Wolstencroft and Smith 79] have suggested the representation 

P(A)/P,nax = 2'^ (A/A 

max ~r -^max 

When K. = 2.25 this curve lies within 1 % of the Serkowski curve with K = 1.15 in the 
wavelength interval 0.22 — 1.40 /im. 
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The values of Pmax in the diffuse interstellar medium usually do not exceed 
10%. The ratios Pj^in:^/ E{B — V) and Pmax/^v determine the polarizing 
efficiency of the interstellar medium in a selected direction. There exist 
empirically found upper limits on these ratios (tsI 

P^^J E{B ~V) <9%/ma.g and P^axMv ^ 3%/mag . (11) 

The mean value of Amax is 0.55 fim although there are directions for which 
Amax is smaller than 0.4/im or larger than 0.8 fxm j8o| (see also Fig. |6]). 
Using observations of about 50 southern stars Whittet and van Breda 



81| established a relation between the parameters of the extinction and po- 
larization curves Ry = (5.6 ± 0.3) Amax, where A^ax is in fim. However, 
further investigations of separate clouds questioned this correlation (e.g., 



67|,|82|,l83|). 

The connection between the coefficient K and the width of the normalized 
curve of interstellar linear polarization W is given by the relation 

W = exp[(ln2/ir)^/2] _ gxp [-(In 2/7^)^/2] _ 
Treating as a third free parameter of the Serkowski curve, Whittet et 



al. j8J| evaluated the dependence between K and Amax on the basis of obser- 



vations for 109 stars 

K = (1.66 ± 0.09)Amax + (0.01 ± 0.05). (12) 

The coefficients of the linear function ( IT2|) for different regions may strongly 
deviate from the average values (see Fig. [6] where the data for the Taurus 
dark cloud and the p Oph cloud are plotted). 

In parallel with the positive correlation between K and Amax, the negative 
correlation between the polarization efficiency Pmax/^v and Amaxfor stars in 

(see also 



separate interstellar clouds and associations is observed [82|, l88l . 
discussion in Sect. 13. 3p . 

The IR continuum polarization for A > 2.5 fim cannot be represented by 
the Serkowski curve with three parameters. The polarization seems to have a 



Another approximation has been proposed by Efimov [80| 

P(A)/P„,ax = [(A/A„,ax) exp(l - A/A 
where the index f3 is proportional to K. 
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Figure 6: The coefBcients K of the Serkowski curve pi))) in dependence on the wavelength 
of the maximum polarization Amax- The data were taken from 80, 82] for the Taurus dark 
cloud and from [85| for p Oph cloud. The HD numbers of stars are marked. Taurus cloud 
1 contains 14 stars with similar positional angles of polarization 0gai. = 145° — 175° (see 
[sij for details). HD 150193 is Herbig Ae/Be star with intrinsic polarization 87|. The 
linear fits for p Oph cloud are shown for cases with and without HD 150193. 



common, universal functional form independent of the value of A^ax and its 
wavelength dependence is given by a power law P{\) oc A"^^'^"^'*^-* 10|, |85 



The UV polarization for 28 lines of sight in the Galaxy has been analyzed 
by Martin et al. 90|] and fitted by a Serkowski-like curve. 

As interstellar extinction and interstellar polarization have different wave- 
length dependencies, the polarizing efficiency P{X)/A{X) has a maximum in 
the near-IR. Note that the polarizing efficiency generally increases with wave- 
length for A ^ 1 /xm. It may be approximated by the power-law dependence 
P/A oc A^. For stars presented in Fig. [7] the values of e vary from 1.41 for 
HD 197770 to 2.06 for HD 99264. 

Variations of the polarizing efficiency in cold dark clouds and star-forming 
regions are of special interest. It was found that in several dark clouds the 
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Figure 7: The polarizing efBciency of the interstellar medium in the direction of five stars. 
Observational data were taken from (extinction) and [9l|, [o^l (polarization). Adapted 
from ^66] . 



rise of polarization with growing extinction was stopped at some value of 
Ay josi 94 1 and that the polarizing efficiency Pma^JAy or Pk/Ak declines 



rapidly with increasing extinction |82l . l83l . l95l l96l . |97[. This fact is usually 



considered as the evidence for the lower efficiency of grain alignment in dark 
clouds in comparison with diffuse clouds [98|. But it is possible to propound 



many other factors influencing t he p olarization degree. They are discussed 
qualitatively by Goodman et al. j93| (see their Table 4). Netherveless, there 
exists the direct evidence for grain alignment in cold dense environments. 



For example, Hough et al. j99|, Il00| . using spectropolarimetry of the 3.1 /im 
and 4.67 /im solid H2O and CO features along the line of sight to Elias 16, a 
fleld star background to the Taurus dark cloud, found that the features were 
polarized. This indicates the presence of multi-layered nonspherical grains 
in molecular clouds because solid CO survive at T < 20 K and solid H2O at 
higher temperatures. 

An error often made is ignoring the foreground polarization. Thi s can 
be well illustrated by the observational data of Messenger et al. 101| who 



analysed the interstellar polarization in the Taurus dark cloud. The authors 
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applied a two-component model and excluded the foreground polarization 
for star HD 29647. Using the data from Tables 1 and 3 of Messenger et al. 



Table 3: Polarization efficiency for stars in Taurus dark cloud*. 



Star 


^max) % 


Av 




Comment 


HD 283812 


6.30 


1.91 


3.30 


cloud 1 (foreground) 


HD 29647 


2.30 


3.32 


0.69 


clouds 1 + 2 


HD 29647 


6.17 


1.41 


4.38 


only cloud 2 (background) 



data from Messenger et al. 10l| 



101| it is possible to estimate the polarizing efficiency in the foreground and 
background clouds. As follows from Table |3l the polarization efficiency in 
cloud 2 (background) is very high (cf. (Eq. ITTl)). Therefore, interpretation of 
the observed polarization instead of the true polarization in the background 
cloud is a large mistake for HD 29647. 

3.2. Interpretation: particles and alignment 

The interpretation of polarimetric observations includes computations of 
the polarization cross sections and their averaging over given particles size 
and orientation distributions. The linear polarization of non-polarized stellar 
radiation passing through a cloud with a homogeneous magnetic field and 
rotating particles can be found as (cf. Eq. (JTj)) 

^(^) = 5Z/ / ^pou (^71, rv, A) n, (ry) dry,,- d/- 100%, (13) 



minj 



Cpoi,(A) = 4r rr kc™-CJX,)m,P)cos2i^dvdojdp, (M) 

where nj{rv) is the size distribution of non-spherical dust grains of the type 
j, ry radius of equivolume sphere (for infinite circular cylinders the parti- 
cle radius r^yx is used). The superscripts TM and TE denote two cases of 
orientatio n of the electric vector of the incident radiation relative to the par- 



ticle axis |102| . The average polarization cross sections Cpoi,j depend on the 



alignment function /(^,/3) with the alignment parameter ^. Here, (3 is the 
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precession-cone angle for the angular momentum J which spins around the 
direction of the magnetic field Lp the spin angle, u the precession angle 
(see Fig. [8]). In general, the particles are assumed to be partially aligned: the 
major axis of the particle rotates in the spinning plane which is perpendic- 
ular to the angular momentum which spins (processes) around the direction 
of the magnetic field. The angle between the line of sight and the magnetic 
field is (0° < < 90°). 



X 




k <f 



Figure 8: Geometrical configuration of a spinning and wobbling prolate spheroidal grain. 
The major (symmetry) axis of the particle O1O2 is situated in the spinning plane NO1O2 
which is perpendicular to the angular momentum J. The direction of light propagation 
k is parallel to the Z-axis and makes the angle a with the partic le sy mmetry axis. The 
angle between the line of sight and the magnetic field is 51. After [103| . 

The cross sections C™'™ = nvy Q™'™ in Eq. f lT^ can be calculated by 
using the light scattering theory for non-spherical particles. To facilitate the 
calculations, particles of simple shapes are usually considered (see Table H]). 
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Table 4: Models used for interpretation of interstellar polar- 
ization. 



Author (s) (year) reference 


Dust grains 


Alignment 






(angle) 


Wilson (1960) [104] 


infinite cylinders* 


PDG 






= 30°, 90°) 


Greenberg et al. (1963-8) 


infinite cylinders. 


PF, PDG 


[39. 105. 106] 


homogeneous spheroids** 




Rogers and Martin (1979) 


homogeneous spheroids 


PF 


[107] 




{a = 90°) 


Hong and Greenberg (1980) 


silicate core— ice mantle 


IDG, PDG 


[108] 


infinite cylinders 




Mathis (1979, 1986) 


infinite cylinders*** 


PDG+ 


[109, no] 


(spheroids) 


(n = 90°) 


Onaka (1980) [111] 


core-mantle spheroids 


PF 

(a = 90°) 


Vaidya et al. (1984) [112] 


homogeneous 


V / 

PF 




spheroids**** 


(a = 45°, 90°) 


Voshchinnikov et al. 


silicate core-ice mantle 


IDG, PDG 


(1986, 1989) [113, 114] 


infinite cylinders 




Mishchenko (1991) [115] 


homogeneous spheroids 


IDG, PDG 


Wolff et al. (1993) [116] 


homogeneous, coated 


PDG+ 




and composite infinite 


{n = 90°) 




cylinders, homogeneous 






spheroids 




Kim and Martin (1994) [117] 


infinite cylinders 


PDG 






(fi = 90°) 


Kim and Martin (1995) [118] 


homogeneous spheroids 


PF, PDG 






(a, n = 90°) 


Matsumura and Seki (1996) 


homogeneous spheroids 


PF 


\m 


and ellipsoids 


{a = 90°) 


Li and Greenberg (1997) 


silicate core-organic 


PDG 


[120] 


mantle finite cylinders 


(fi = 90°) 


Vaidya et al. (2007) [121] 


silicate spheroids with 


PF (a = 45°, 




graphite inclusions 


60°, 90°) 


Wurm and Schnaiter (2002) 


dust aggregates 


PF 
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Table 4: (Continued.) 



Author (s) (year) reference Dust grains 



Alignment 
(angle) 



122 



Voshchinnikov et al. (1990- 
2010) [66. 123. 124. 103] 
Draine and Fraisse (2009) 
[1251 



consisting of 4-64 
monomers 

homogeneous spheroids 
oblate spheroids 



(a = 90°) 

IDG, PDG 

PDG++ 
(fi = 90°) 



PF — picket fence orientation 

PDG — perfect Davis-Greenstein (2D) orientation 

IDG — imperfect Davis-Greenstein orientation 

a, (fi) — angle between the line of sight and the direction of grain alignment in 
the case of PF (PDG) orientation 
* in Rayleigh-Gans approximation; 

** in Rayleigh approximation; 

*** efficiency factors tabulated by Wickramasinghe 126l | were taken; polarization 
for prolate spheroids was com puted as for cylinders; 
**** the figures of Asano [l27| were used; 

+ large silicate grains are assumed to be perfectly aligned (see Eq. (fTSl) ): 

only grains with sizes r > rcut are assumed to be perfectly aligned (see Eq. (I19p ). 



Ear l y models dealt with homogeneous infinitely long circular cylinders [105 



Toeiliogl. Eol pl. This is the simplest model of non-spherical particles. So 
lution to the light scattering problem for infinite cylinders was obtain ed by 



the separation of variables method in the cylindrical coordinate system 128 



Later, more advanced models wit h sil icate core-ice mar itle c ylindrical parti- 
cles based on the solution from 129| were developed lOSL Ill3l . Ill4j . The 
progress in the light scattering theory allowed one to apply the model of 
homogeneous prolate and oblate spheroids of different size and shape for cal- 
culation^ o£_thep^^ efficiency, visual and UV polarization (see, e.g.. 

Spheroidal particles are characterized by the aspect 
ratio a/h where a and h are the major and minor semiaxes. The optical 



* Another simplification was the use of the wavelength-independent refractive index. 
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properties of spheroids can be found by using different teclinique. Tlie most 
popular methods wid ely apphe d in astronomical model ling a re the separation 



of variables method |l24j. 
dipole approximation 131 



30[, the T- matrix method |115| and the discrete 
Comparison of the methods and benchmark 



results are given in 132i 133| pl. More complicated non-spherical particles 
(coated spheroids, ellipsoidal parti cles, composite spheroids) were considered 



so far for illustrative calculations 111. 119. 121 



The polarization cross sections must be averaged over rotations taking 
into account an alignment mechanism. According to standard concepts (l|, 
13J], the alignment of interstellar grains may be magnetic or radiative. A very 
popular alignment mec hani sm is the magnetic alignment (Davis-Greenstein 
(DG) type orientation 135| ) based on the paramagnetic relaxation of grain 
material containing about one percent of iron impurities. For the imperfect 
Davis-Gre enstein (I DG) orientation, the alignment function /(^,/3) can be 



written as 



ioi.[Ti^ 



/(e,/3) 



^ sin (3 



ie cos2/3 + sin2/3)3/2- 



(15) 



The parameter ^ depends on the particle size ry, the imaginary part of the 
grain magnetic susceptibility x" (= >^^^d/Td, where Ud is the angular velocity 
of grain), gas density Ug, the strength of magnetic field B and dust (Td) and 
gas (Tg) temperatures 



ry + So{Td/Tg 
rv + ^0 



where = 8-23 10^^ 



nJT' Td 



yUm. 



(16) 



If the grains are not aligned ^ = 1 and /(^,/3) = sin/3; in the case of the 
perfect rotational orientation ^ = 0. Unfortunately, only a limited number of 
models including the combined particle size/shape/orientation analysis have 
been developed (see Table H]). 

For simplicity, many investigators assumed that the direction of the mag- 
netic field (direction of gra in al i gnnaent) was perpendicular to the line of 
sight, i.e., a, n = 90° (e.g., fTT0l.flT8l.E25l I Frequently, non-rotating parti- 
cles of the same orientation are considered (see Table H]). In this case, called 
the "picket fence" (PF) orientation, there are no integrals over the angles 



^see also Database of Optical Properties (POP) : [http : //www . astro . spbu . ru/DOP 
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(jj and /3 in Eq. ( IT^ . The polarization degree is proportional to the polar- 
ization cross-section P oc Cpoi = l/2[C™(f2) — C™(f2)], where Vt = a and 
f{^,(3) = 5{a^. The dichroic polarization efficiency is defined by the ratio 
of the polarization cross-section (factor) to the extinction one 



pol 



a 



TM 
ext 



■ 100% 



Q - Q 



TE 

■ 100%. 



(17) 



PF '-"ext '-"ext "-"ext Vext Vext 

A more complicated case is the perfect rotational (2D) orientation (or perfect 
Davis-Greenstein orientation, PDG) when the major axis of a non-spherical 
particle always lies in the same plane. For the 2D orientation, integration is 
performed over the spin angle ip only and /(C,/3) = S{(3). 

Polarization produced by perfectly aligned particles is much larger than 
that observed (cf. Figs. |9] and [TO] with Eq. f lTTj) ). Nevertheless, the models 
with the PF or PDG orientation are useful for investigations of the normal- 
ized polarization (the Serkowski curves) as the wavelength dependence of 
polarization is only slightly infiuenced by the particle refractive index, size 
or shape (cf. left panels in Figs. |9] and [TOj) . As a crude approximation of the 
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Figure 9: Wavelength dependence of the polarization efhciency for homogeneous spheroids 
consisting of astrosil for the PF orientation (see Eq. ([T7|) ). The effect of variations of the 
particle size (left panel), type and orientation (right panel) is illustrated. Adapted from 



polarization efficiency for part icles with the IDG orientation, the following 
relation is used (see, e.g.. 



7^ sin^ Q 



IDG 



PF 



^S{z) is the Dirac delta function. 
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Figure 10: Wavelength dependence of the polarization efhciency for homogeneous rotating 
spheroidal particles of the astronomical silicate. The effect of variations of particle size (left 
panel), and degree of alignment (right panel) is illustrated. The open circles and squares 
show the observational data for stars HD 24263 and HD 99264, respectively. Adapted 
from [66| . 



where IZ = l/2(3(cos^ (3) — 1) is the Rayleigh reduction factor [39|, Il36l | and 
denotes the ensemble average. It should be emphasized that application 
of this approximation as well as the Rayleigh reduction factor can lead to 
misinterpretation of observational curves -P(A). 

In the case of the IDG mechanism, smaller grains are aligned better than 
larger grains (see Eq. ( IT6|) ). However, the models with an opposite ty pe of 



orientation of small and large particles have been also suggested. Mathis |110 
assumed that rotating silicate grains were perfectly aligned if they contain 
at least one super-paramagnetic inclusion. Carbonaceous grains and silicate 
grains without inclusions are randomly oriented in space (3D orientation). 
The probability of perfect alignment is 

f{rv,r'y) = l-exp{-rv/r'y)\ (18) 

Draine and Fraisse jl25 | considered the model of silicate and amorphous 



carbon spheroids with randomly oriented small particles and perfectly aligned 
large particles. In this case the alignment function is size dependent 

for ry>ry,,^,, ^^^^ 

where ry^cut is a cut-off parameter. 

Computations made by Das et al. 103l | (see their Fig. 2) demonstrate 
that the observational data can be fitted by using the models with different 
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alignment functions (e.g., given by Eqs. ( IT^ . (fTHj) or (fT^ ). especially if a 
more complex size distribution function as discussed in Sect. 12.21 is chosen. 
However, this complicates the model. To avoid the models with many pa- 
rameters new ideas about the nature of polarizing grains and physics of grain 
alignment should be included. 

The DG mechanism of the paramagnetic relaxation requires a st ronger 
magnetic field than average Galactic magnetic field (~3 - 5 /iG; |l37 |). Be- 
cause of this problem, it has been suggested that the polarizing grains contain 
small clusters of iron, iron sulfid es, or iron oxides with super-paramagnetic 
or ferromagnetic properties 138l | . This leads to an enhancement of the imag- 
inary part of the magnetic susceptibility of grain material x" by a factor 10 
- 100 and alignment can occur throug h the DG mechanism. This scenario is 
supported by laboratory experiments |l39 . 140 |. A significant enhancement 



of x" is also possible in mixed MgO/FeO/SiO grains |141] or in H2O ice 
mantle grains containing magnetite (Fe304) precipitates [ 142[ Il43 . 



Another possibility to align interstellar grains is the radiative torque 
alignment (RAT alignment). It arises from an azimuthal asymmetry of the 
light scattering by n on-spherical particles. Magnetic inclusions can en hanc e 



RAT alignment |144| . The theory of RAT alignment is well developed [145 



Rece nt observations of interstellar polarization in the vicinity of luminous 



stars |146l . Il47l . Il48| have been used for confirmation of the RAT alignment 
mechanism. However, the discussed models are phenomenological, they are 
not based on correct light scattering calculations of interstellar polarization. 
One of the reasons is that the alignment function for the RAT mechanism 
is unknown. Another reason is a requirement of advanced light scattering 
methods bec ause fast rotation can only occur for grains of very specific (he- 



lical) shape 145l . Il49| . This is highly improbable from the point of view of 
grain growth in the interstellar medium. 

Since both magnetic alignment and radiative alignment depend on iron 
inclusions, we can expect that polarization and/or polarization efficiency 
should increase with the growth of iron fr actio n in dust grains. This idea 



was investigated by Voshchinnikov et al. |150l | by using available data on 



interstellar polarization and element abundances previously compiled in [29 



It was suggested that the interstellar polarization was probably related to the 
amount of iron in dust grains. Assuming that all silicon and all magnesium 
are embedded into amorphous silicates of olivine composition (Mg2xFe2-2a;Si04, 
where x = [Mg/H](i/(2[Si/H]d) as is a part of iron. The remaining part of Fe 



25 



can be found as 

[Fe(rest)/H], = [Fe/H], - (2 [Si/H], - [Mg/H],). (20) 
As indicated in Fig. [TT] (left panel), there is a negative correlation between 
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Figure 11: Interstellar polarization in dependence on remaining dust phase abundance of 
Fe as given by Eq. (PU]) (left panel) and dust phase abundance of silicon (right panel). 
Halo stars with |6| > 30° and disk stars with |6| < 30° and low {E{B - V) < 0.2) and 
high {E{B — V^) > 0.2) reddening are shown with different symbols. Number of stars is 
indicated in parentheses in the legend. Adapted from 150j . 

the polarization degree P and the amount of remaining iron. This is incon- 
sistent with the common suggestion about the great role of iron-rich grains in 
the production of polarization. Since P is proportional to the column density 
of polarizing grains, we can conclude that the increase of the iron content in 
non-silicate grains does not enhance polarization. 

Because in calculating [Fe(rest)/H]^ we removed all Si and Mg and a 
part of Fe from the dust phase, we expect a positive correlation between 
the polarization and the abundances of the eliminated elements. Ther e is 
only a weak correlation between P and [Fe/H]^ or [Mg/H]^ (see jl50| for 
more discussion) and a strong correlation between P and [Si/H]^ (Fig. [HI 
right panel). Therefore, it can be established that polarization is more likely 
produced by sil icates. These findings are evidence in favour of the assumption 
of Mathis that only the silicate grains are aligned and contribute to the 
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observed polarization, while the carbonaceous grains are either spherical or 
randomly aligned. Another verification of this suggestion is the absence of 
any correlation between the polarization effic iency P/ E[B — V) or P/Ay and 
dust phase abundances of elements (see |l50[ |). This is because dust grains of 
all types (silicate, carbonaceous, iron-rich, etc.) contribute to the observed 
extinction, while only the silicates seems to be responsible for the observed 
polarization. Thus, the absence of correlation between Ry and Amax (see 
discussion in Sect. 13. ip can be easily understood. These discoveries can be 
explained if the silicate grains aligned by the radiative mechanism are mainly 
responsible for the observed interstellar linear polarization^. 

Analysing models presented in Table HI it is possible to say that the major 
part of models includes perfect grain alignment and one angle of alignment 
and merely a few of them with the IDG orientation can be used for inter- 
pretation of observations of individual stars. Indeed, the previous modelling 
of interstellar polarization was mainly focused on the explanation of the av- 



erage wavelength dependence (Serkowski curve) |l09 . lid . 117 . 118]. Only 



Li and Greenberg [69| applied their model of coated cylinders to explain 
the no rmal ized polarization curve in the direction of HD 210121 and Das 
et al. 103l | interpreted interstellar extinction and polarization observations 



of seven stars using a mixture of carbonaceous and silicate spheroids. This 
fact causes deep dissatisfaction because a great amount of observations of 
interstellar polarization in different areas exists and continuously grows. 

3. 3. Interpretation: dust grains and magnetic field in the Taurus dark cloud 

In this section we present the quantitative interpretation of observations 
of interstellar polarization for a gro up o f star^. Our model of spheroidal 



grains with imperfect alignment [66|, Il03l | is supplemented by the subroutine 
calculating three parameters of the Serkowski curve. We also assume that 
polarization is mainly produced by silicate grains and the degree of alignment 
of carbonaceous grains is small (see discussion in previous section). 

As an example we refer to the Taurus dark cloud (TDC) — the complex of 
interstellar clouds where active star formation is in progress. This complex 



^Polarization in IR features also supports the idea of separate populations of polarizing 
(silicate) and non-polarizing (carbonaceous) grains. This follows from the observed polar- 
ization of silicate features at 10 Atn i and 18 /zm and the lack of polarization in the 3.4 /zm 
hydrocarbon feature (see [l5ll Il52| and references therein) . 

^More detailed discussion can be found in [si'l. 
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lies sufficiently far from the Galactic plane (6 ~ —15 
^ 1 40 pc and comprise several dozens dark nebulae, clouds and clu mps 



at a distanc e of 



154| . It suffers negligible foreground and background extinct ion 1 155 



and is th e subject of numerous investigations of molecular gas |1571 1158 



15 



15 



160, 161 



The polarimetric observations of several hundreds stars in the region 
around the Taur us dark c loud complex have been performe d at visual and 
red wavelengths 160l . Il62l . Il63j and in the H and K bands 95|, Il64j . Il65l Il66 



The obtained polarization maps give an insight into the structure of the 
plane-of-the-sky component of interstellar magnetic field but are of little 
use for studies of grain properties and alignment. To find these latter, the 
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wavelength dependence of polarization must be involved. Whittet et al. 
presented the polarimetric and photometric observations of 27 stars in a wide 
spectral range in the TDC {I ^ 170°-^ 176°, b ^ -10°^ -17°) and calculated 
the fit parameters of the Serkowski curve Pmax, Amax and K as well as the val- 
ues of Ry- Our initial analysis of these data is base d on the two-component 



model of Messenger et al. |l01| and Whittet et al. |l67l | (see also discussion 
of Table Elin Sect. 13. ip . It made possible to form two groups of stars with 
relatively uniform distribution of positional angles of polarization: cloud 1 
(14 stars, ^gai. = 145° - 175°) and cloud 2 (13 stars, ^gai. = 2° - 40°). 

Firstly, we should concentrate on stars in cloud 1 located inside or be- 
hind the "diffuse-screen" component of the TDC. These sta rs are distributed 
around Heiles Cloud 2 {I ^ 174?4, b ^ -13?4, area 15.8 pc^ [l6l|) — a dense 
condensation containing TMC-1 (Taurus molecular cloud). Observational 
data for stars in the cloud 1 plotted in Figs. [121 ESI indicate positive correla- 
tion between K and Amax (the width of the polarization curve W decreases 
with increasing Amax) and negative correlation between the polarization effi- 
ciency -Pmax My and Amax- The ffist dependence is well known (see Fig. Eland 
discussion in |8J]). Its qualitative explanation is a systematic reduction in 



the relative number of small, aligned grains in regions of high Amax |82l . |90 



The sole quantitative modelling was attempted by Aannestad and Greenberg 
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However, to reduce computational efforts they ruled out the integra- 
tion over the angle u in Eqs. (I14p which led to wrong results for W 114 . 
A systematic trend toward smaller polarizing efficiency for larger Amax was 

It has been 



82, 88 



detected for nearby stars closely located on the sky 
interpreted as a result of the decrease of t he a ngle Q between the direction 
of the magnetic field and the line of sig ht fnl . 



The results of theoretical modelling are shown in Figs. [T21 ESI by open cir- 



28 




Figure 12: The coefficients K of tlie Serkowski curve (fTO|) in dependence on tlie wavelength 
of maximum poiarization Amax for 14 stars in cloud 1 in Taurus with the similar positional 
angles of polarization. The numbers of stars correspond to the increasing HD numbers 
of stars in Fig. H i.e. 1 = HD 28225, . . . , 14 = HD 283879. The Pearson correlation 
coefficient between K and Amax is given. Open circles with line show model calculations 
with different alignment angles VI = 15°(15°)90°. The values of VL increase from right to 
left as marked for the top model. Left panel illustrates the effect of variations of the power 
index q in the power-law size distribution {q varies from -0.5 to -3). Right panel illustrates 
the effect of variations of the lower cut-off rv,min in the power-law size distribution (ry^min 
varies from 0.03 /xm to 0.12 ^va). Other model parameters are: prolate spheroids, a/b — 3, 
''v,max = 0.35 /im, rv,min = 0.07/^m (left panel), q = —2 (right panel). 



cles connected with solid line. The particle shape was fixed: prolate spheroids 
with the aspect ratio a/b = 3. Under assumptions made the main parameters 
of our model infiuencing the polarization are: the lower and upper cut-offs 
rv,mm and ry,max and the power index q in the power-law size distribution for 
silicate particles and the degree (parameter S^^, see Eq. (fT6!) ) and direction 
(angle Q) of alignment. Left and right panels in Fig. [12] illustrate the vari- 
ations of the index q and lower cut-off rv,ram, respectively. The rise of these 
parameters may be associated with the growth of dust grains by coagulation 
(g) or accretion {rv^mm)- In both cases the mean size of grains is bigger at 
the right upper corner of Fig. [12] in comparison with its left bottom corner. 
It is interesting that the stars NN 1, 5, 6, 8, 11, 14 located at the right upper 
corner apparently are embedded in Heiles Cloud 2 or situated at its boundary 
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Figure 13: The polarizing efficiency -Pmax/^y in dependence on the wavelength of maxi- 
mum polarization Amax for 14 stars in cloud 1 in Taurus. The Pearson correlation coeffi- 
cient between Pmax/^y and Amax is given. Open circles with line show model calculations 
with different alignment angles fl = 15°(15°)90°. The values of ft increase from bottom 
to top as marked for the right model. Left panel illustrates the effect of variations of 
the power index q in the power-law size distribution {q varies from -0.5 to -3). Right 
panel illustrates the effect of variations of the lower cut-off ry^min in the power-law size 
distribution {rv,min varies from 0.03 fiui to 0.12 fim). Other model parameters are: prolate 
spheroids, a/b = 3, rv.ma.x = 0.35 /im, r^/^min — 0.07 /im (left panel), q = —2 (right panel). 



(see Fig. [T^ . Note also that variations of alignment parameters and Q 
only do not allow explaining the observed correlation between K and Amax- 

The opposite situation occurs with the correlation between the polariza- 
tion efficiency PmaxMv and Amax (Fig. [I3]). In this case to explain both 
high and low values of Pmax/^v it is not sufficient to change the grain size 
only, variations of parameters 6^"^ and Q must be taken into account. The 
theoretical points plotted in Fig. [13] were obtained for = 3 fim. Smaller 
values of do not reproduce the data for directions with high polarization 
efficiency (stars NN 2, 3, 10). It is evident that the obtained grain parameters 
are model-dependent (e.g., it is possible to use particles of another shape). 
However, with our model of interstellar dust the trends in variations of the 
grain size, shape, alignment can be determined. 

We attribute the variations of Pmax/^v in Fig. [T3]to the changes in the 
direction of magnetic field (direction of grain alignment). Stars can be di- 
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Figure 14: Linear polarization map of the cloud 1 in Taurus containing 14 stars with similar 
positional angles. The lengths of the lines are proportional to the percent polarization. 
The sizes of the circles are proportional to the alignment angle fi. The dashed contour 
represents the regions with different vi sual extin ction: 1™6 < Ay < 19™6 inside contour 
and 0™4 < Ay < 1™6 outside con tour 156 . 169l |. The cross corresponds to the clump in 



TMC-1 where Turner and Heiles |170| have searched the C4H Zeeman effect. The open 
circles show the positions of dark clouds where the OH Zeeman effect has been observed 

ml- 



vided into three groups with close values of Q (see right panel in Fig. [T3|) . 
These groups of stars are shown in Fig. [14] by filled circles of different sizes. 
This Figure gives the positions and polarization of stars in cloud 1 and ap- 
proximate contours of Heiles Cloud 2. It is intriguing that all six stars with 
the larger values of Q (the magnetic field is mainly perpendicular to the line 
of sight) are closely located on the sky outside of Heiles Cloud 2. Other eight 
stars where the magnetic field is significantly tilted to the line of sight are 
situated near the boundary of the cloud. An indirect support t o our interpre- 
tation is provided by Zeeman observations. Turner and Heiles 170l | obtained 
an upper limit for magnetic field toward the cold dense TMC-1 cyanopolyyne 
peak core Bii = B cosVt = 14.5 ± 14 fiG. A possible reason for this result is 
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tha t the magnetic field is directed close to the plane of the sky. Crutcher et 



al. |171| compiled the Zeeman data for many dark clouds. Several of them 



located near or inside Heiles Cloud 2 are marked as large open circles in 
Fig. [TH In almost all cases the data are within la error. Only for B217-2 
(close to the positions of stars NN 5, 6) the component of the magnetic field 
parallel to the line of sight has been detected {B\\ = 13.5 ± 3.7yuG). 

The structure and evo lutio n of magnetic field in H eiles Cloud 2 have been 



discussed by Heyer et al. [1601 ] and Tamura et al. jl64| . It was suggested that 



the contraction and formation of the cloud from the placental Taurus dark 
cloud with homogeneous magnetic structure occurred along the magnetic 
field. The gas motions in Heiles Cloud 2 can be described in terms of a 
rotating ring with the rotation axis coinciding with the magnetic field 157 . 
Our findings infer that the magnetic field is tilted with respect to the line of 
sight at the boundaries of Heiles Cloud 2, i.e., magnetic field has a spindel- 
like structure. Apparently, the area disturbed by contraction of Heiles Cloud 
2 terminates at the place where the magnetic field is almost perpendicular 
to the line of sight (stars NN 2, 3, 9, 10, 12, 13). 



4. Conclusions 

All modern investigations of cosmic dust in protoplanetary disks, dense 
clouds, distant galaxies are based on the modelling of "classical" observations 
of the interstellar extinction and polarization. Our examination shows that 
the interpretation of these basic observations still remains incomplete. New 
generation models should include a consideration of interstellar abundances 
in given directions and accurate treatment of light scattering by physically 
feasible particles with realistic alignment. 
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